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ABSTRACT 

The  fundamental  equations  governing  the  process  of  isotope 
separation  by  thermal  diffusion  are  derived.  The  equations  are  applied 
to  the  behavior  of  a  single  column  in  the  cases  of  equilibrium, 
stationary  flow  and  approach  to  equilibrium. 


UNCLASSIFIED 


INTRODUCTION 

1.  The  s  eparafcion  of  isotopes  by  the  method  of  thermal  diffusion 
has  assumed  a  practical  significance  in  connection  with  the  problem  of 
obtaining  concentrated  uranium  235,  an  element  which  is  capable  of 
releasing  large  amounts  of  energy.  .hile  the  work  on  this  problem  is 
primarily  of  an  experimental  nature,  theoretical  work  can  play  an 
important  role  in  guiding  the  work  of  the-  experimenter  and  in  correlating 
his  results. 


2,  The  present  paper  deals  with  the  theory  of  the  thermal  diffusion 
method.  In  the  published  literature  there  are  a  number  of  papers  on 
the  theoretical  treatment  of  isotope  separation  by  thermal  diffusion. 


among  which  one  might  single  out  for  mention  those  of  Furry.  Jones  and 
Onsager^(to  be  referred  to  as  F.  J.  0.),  Haldmann^ ,  Bardeen  ,  and  Debye^. 
However,  all  of  these  but  the  last  are  concerned  with  gases  rather  than 
liquids.  The  paper  of  Debye  is  not  Satisfactory  because  of  the  rough 


approximations  introduced,  while  the  papers  dealing  with  gases  make 
some  assumptions  which  are  probably  not  valid  for  liquids,  (incidentally 


it  is  interesting  tc  note  that  since  1940  no  papers  on  isotope  separation 


by  thermal  diffusion  have  appeared  in  German  publications,  although 
there  had  been  numerous  papers  prex^iously  .) 


3.  Because  there  exists  no  satisfactory  theoretical  treatment  of 
isotope  separation  by  thermal  diffusion  in  liquids,  it  seems  desirable 
to  attempt  such  a  treatment.  The  purpose  of  the  present  paper  is  to 
discuss  the  case  of  a  single  liquid  diffusion  column,  introducing 
assumptions  and  approximations  which  are  appropriate  to  the  type  of 
column  being  used,  and  to  obtain  results  which  can  be  directly  applied 
to  it. 


APPARATUS  AND  PROCESS 

4.  In  order  to  be  able  to  treat  the  problem  theoretically,  we 
must  first  have  a  clear  picture  before  us  of  the  system  to  be  investi¬ 
gated.  Some  idealization  will  be  necessary,  to  be  sure,  but  it  is 
important  to  try  to  retain  the  essential  features. 

5.  We  shall  consider  a  column  to  consist  of  two  vertical 
concentric  tubes,  spaced  a  distance  a  apart,  each  of  length  L,  between 
which  the  diffusing  liquid  is  contained,  in  the  annular  space  of  mean 
circumference  b.  The  surfaces  of  the  tubes  in  contact  with  the  liquid 
will  be  referred  to  as  the  "walls”.  One  wall  is  maintained  at  a  low 
temperature  Tl,  the  other  at  a  high  temperature  T2*  The  mass  of  the 
fluid  contained  between  the  walls  will  be  denoted  by  lic.  At  the  top 
of  the  column  there  is  a  reservoir  filled  with  liquid  having  a  masa’ 
Mp.  At  the  bottom  of  the  column  there  is  an  outlet  which  iray  be  con¬ 
nected  to  another'  reservoir,  or  else  kept  closed, 

6.  Because  of  the  temperature  difference  between  the 'walls,  a 
con  ■action  current  sets  in,  the  liquid  near  the  hot  wall,  having  the 


>■ 

oo 
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lower j  density,  flowing  upward,  while  that  near  the  cold  wall  flows  down. 

At  the  same  t ime  thermal  diffusion 'takes  place,  the  molecules  of  the 
lighter  isotope  diffusing  from  the  cold  toward  the  hot  wall.  The  result 
is  that  there  is  an  exchange  of  molecules  between  the  Upward  convection 
stream  and  the  downward  one,  by  which  the  upward  stream  acquires  an 
excess  of  the  lighter  isotope  and  carries  it  to  the  top  of  the  column. 

The  concentration  of  the  lighter  isotope  at  the  top  of  the  column,  or 
in  the  reservoir,  gradually  increases,  while  that  at  the  bottom  decreases, 
unless  the  bottom  is  connected  to  a  large  reservoir  or  to  another 
column,  so  that  a  constant  concentration  is  maintained* 

7*  If  this  process  is  allowed  to  continue,  a  condition  of 
equilibrium  will  be  approached,  at  which  there  will  be  a  definite 
relation  between  the  concentrations  at  top  and  bottom.  It  is  possible 
to  stop  the  process  at  some  stage,  remove  the  enriched  liquid  from  the 
reservoir  and  start  the  process  over  again.  Such  a  procedure  is  referred 
to  as  intermittent  operation.  On  the  other  hand,  it  is  possible, 
beginning  at  a  certain  stage,  to  draw  off  a  continuous  stream  of  enriched 
liquid  from  the  top  of  the  column  •  This  procedure  is  called 
continuous  operation.  Under  working  conditions,  the  rqte  at  which  the 
isotope  is  drawn  off  from  the  top  of  the  column  is  equal  to  that  at 
which  it  is  carried  to  the  top  of  the  column  by  convection,  so  that  one 
has  a  stationary  condition  within  the  column. 

8.  In  the  mathematical  treatment  of  a  column  such  as  has  been 
described,  one  can  take  into  account  the  fact  that  the  distance  between 
the  walls  a  is  very  small  compared  to  the  circumference  b,  by  neglecting 
the  curvature  of  the  walls  and  treating  them  as  plane  surfaces.  In 
discussing  the  convection  and  diffusion,  one  can  simplify  the  calculat¬ 
ions  considerably  by  neglecting  the  small  regions  near  the  ends  of  the 
column  where  the  temperatures  vary  from  the  values  prevailing  elsewhere 
and  where  the  convection  streams  curve  back  to  reverse  their  directions. 
These  regions  can  be  neglected  because  their  dimensions  a  re  generally 
very  small  compared  to  the  usual  length  of  a  column.  (F.J.O.^-) 

9«  Finally,  since  a  column  in  a  practical  case  will  produce  a 
relatively  small  change  in  the  isotope  concentration,  one  can  represent 
the  concentration  very  often  as  a  power  series  in  the  distance  up  from 
the  bottom  of  the  column.  The  series  will  usually  converge  so 
rapidly  that  only  a  small  number  of  terms  are  required  for  reasonable 
accuracy.  Such  a  power  series  expansion  frequently  reduces  the  work 
of  mathematical  calculation  considerably. 

C  OBJECT  I  ON 

10.  Since  we  are  dealing  with  isotopes  having  a  small  percentage 
difference  in  atcmic  weight,  the  two  kinds  of  molecules  being  separated 
are  very  similar  in  their  properties,  and  the  presence  of  thermal 
diffusion  does  not  alter  the  behavior  of  the  liquid  to  any  marked 
extent*  Hence  one  can  first  treat  the  convection  of  the  liquid  in  the 
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.absence  of  diffusion  and  afterwards  treat  the  diffusion  as  influenced 
by  the  convection* 

11.  Let  us  suppose  that  we  have  the  liquid  between  two  parallel 
vertical  walls  a  distance  aap.art,  of  breadth  b  nnd  height  L.  Let  us 
take  a  coordinate  system  with  the  X-axis  at  right  angles  to  the  plates 
and  the  Z-nxis  vertical,  so  that  0 ^  X  $  ft.,  ;>  5  "i  -J-  L  . 

Since  we  shall  consider  the  case  in  which  none. of ^quantities  to  bo 
discussed  varies  in  the  direction  of  the  third  coordinate  axis,  the 
latter  will  not  be  used.  Let  us  take  the  ’’cold"  wall,  as  the  one  at 
x  -  0,  and  let  its  temperature  be  Tq»  while  the  "hot"  wall*  at  x  s  ja, 
has  a  temperature  T2»  so  that  nT  m  Tg  —  Tq>  0.  Let  jp  bo  the  pressure, 
g  the  acceleration  of  gravity,  and  Q  tho  heat  flow  p"er  unit  area  of 
pla  te«  Fur  the  more  let  o  ,  X  and  p  denote  the  density,  thermal  conduce 
tivity,  and  viscosity  of  the  liquid,  as  functions  of  the  temperature 
T»  Debye^  assumed  the  last  three  quantities  to  be  constant,  corres¬ 
ponding  to  a  small  value  of  at*  However,  these  quantities  may  change 
rapidly  with  the  temperature  and  therefore,  for  any  appreciable  value 
of  aT,  the  dependence  on  temperature  should  be  taken  into  consideration. 
On  the  other  hand*  because  of  the  similarity  of  the  two  kinds  of 
molecules,  we  shall  assume  that  p , A  and  p  do  not  do pond  on  their 
cone ontr at ions.  Finally,  from  the  conditions  of  the  problem,  we  can 
take  these  quantities  independent  of  Z# 

12.  To  calculate  the  convection  flow,  we  follow  F.  J#  0*  in 
detc mining  the  temperature  distribution  on  the  basis  of  the  boat 
conduction  alone*  Except  near  the  edges  of  the  plate,  one  can  write. 


<4  f 

a  = 


(i) 


or,  on  integrating. 


(2) 


where 


1 


X  = 


A  J  T. 


(3) 


13*  F.J.O.  showed  that  tho  convection  flow  can  be  treated  as 
lamellar*  To  a  sufficient  accuracy,  one  can  use  the  hydrodyamical 
equation  for  steady  viscous  flow  in  tho  form 


(V-  /-<■  v)  v  -  Vf  -fi  ,  (4) 

where  v,  the  velocity  vector,  is  taken  in  the  z-direction,  but  inde¬ 
pendent  of  z,  while  g,  the  gravitational  acceleration  vector  is  taken 
in  tho  (-z)"**direction  and  of  magnitude  g.  Eq.  (4)  gives  for  the 
x-oemponentj 


is: 

<y> 


m 

C2 
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(4) 


:2L 

Tx 


and  for  the  z  component j 

^  1  =  U  +  >  6‘  (5) 

7)X  ^flX  / 

Since^  in  (5)»  all  terms  ether  than^p  depend  only  on  x,  we  see  that 

=  B,  (6) 

r* 

where  3  is  a  constant,  and  the  equation  can  be  written  by  means  of  (1) 
as  ” 

Q2  d  (  fx  d?  \  =  B  -*  rg-  (7) 

T^{  X  dv> 

v  •  / 

This  equation  is  to  be  solved  subject  to  the  boundary  conditions 

V(Tq)  s  V(T2)  =  0  (8) 

The  solution  of  (7)  can  be  obtained  by  quandrfctures  and  will  depend 
linearly  upon  the  constant  B»  The  value  of  this  constant  will  depend 
on  the  total  flow  of  liquidT  It  will  be  seen  from  (7)  that  for  given 
values  of  B,  Tq  and  T2,  v(T)  for  a  fixed  value  of  T  will  vary  as  l/Q^  or 
by  (1),  as  a~. 

DIFFUSION 

14.  Let  the  two  types  of  molecules  in  the  liquid  be  referred 
to  by  subscripts  1  and  2,  so  that,  for  example,  C-^  and  C2  are  the 
relative  particle  concentrations,  (cq-f--C2  =  1),  and  let  1  represent 
the  species  which  it  is  desired  to  concentrate,  i,e«,  the  lighter  one. 
The  equation  of  diffusion  can  be  written  in  the  form 

ci($i  - "f)  =  -  DC?9q  CVT,  (9) 

where  the  first  term  on  the  right*»hand  side  represents  ordinary  diffuse 
ion,  while  the  second  term  represents  thermal  diffusion.  In  the  case 
of  ordinary  diffusion,  the  diffusion  coefficient  D  (often  denoted  by 
Dq2)  known  to  be  a  function  of  the  temperatureT  It  will  be  assumed 
that  it  is  independent  of  the  concentrations  and  c2,  an  assumption 
which  is  justified  by  the  similarity  of  the  two  species  of  molecules. 

15,  The  coefficient  C  of  the  thermal  diffusion  term  must  depend 
not  only  on  the  temperature’,  but  also  on  the  concentrations.  Since 
there  is  no  thermal  diffusion  if  either  species  alone  is  present,  one 
can  assume  C  to  contain  a  factor  cqCg,  Let  us  write 
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(10) 


C  -  qcxc2  (10) 

D 

on.  CM  expect  that  q  will  depend  oh  T  and  that  it  will  be,  in  general. 

tivo1ape^J^7r)e6vfUn0tl°!'i°f  concentri,tlons •  In  particular,  if  the 

Wie?  tZi Twn!\81"U*r  10  their  it  is  reasonable  t, 

eiieve  that  q  will  be  very  nearly  independent  of  c  and  c  .  This  is 

"  U  th*  ?“•  tj* 

contrary  that  o  in  aS8ume»  until  there  is  evidence  to  the 

this  fS^tion  if  wo  rf  *S  *  function  only  of  the  temperature.  What 
has  been  l  L  to°W  at  /resent*  I*  the  case  of  gases,  it 

this  is  the  function  »  f°ra  **/?»  where  c^is  nearly  a  constant,  and 

argunents^  have  led  to  eithera /^or  J?2  V9ry  r0Ugh  the®retical 

assumed.  It  seems  w*  < orVTS  depending  on  the  model 

leave  open  the  question* of^t  STd^e SE/Sf' ™  * Tis^'  * 
given^  Fr0m  (9)  ^  (10)  the  flux  density  secies  1  is 

1  */,ClVl  s  y^lv  +  ©(-Vcj-^  qcic2 7T)J  ,  (11) 

so  that  we  cm  write 


"bt  --7.^,  (!2) 

such* a?  thet^f .  }* . .  WB  are  dealin6  with  a  stationary  condition, 
member  *  is^s  or  and  ^  "  °f  °P-ation,  the  left-hand 

V.Hi  :  0  (13) 

As  pointed  out  by  F.J.O.1,  the  condition  will  be  essentially 

™se™?T  the  approach  to  equilibrium  provided  the  end- 

reservoirs  are  sufficiently  large. 

STATIONARY  CONDITION 

17.  Substituting  (11)  into  (13),  we  get 

if  /  ^  T\]  ^C,  ^c, 

Tx  LfD(-  57  f  qc^^xjj  r  0,  (14) 

or  if  we  replace  x  by  T,  by  means  of  (1), 


^  i  Z70 

Tf  l/Tr 


( p'  qcici]  =  £  1; v ^  -*&] (i6) 


The  boundary  conditions  at  the  walls  are  given  by 

=  0  (x  z  o,a). 
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18*.  In  as  lving  these  equations,  wo  shall  not  look  for  an  ex¬ 
plicit  solution  giving  the  value  of  c^  for  each  value  of  T  and  z; 
rather  wo  shall  follow  F.J.O*  in  seeking  information  about  tho  average 
value  of  c-^  across  a  section  z-eonstant,  for  various  values  of  z, 
since  this  is  the  quantity  which  is  of  the  most  direct  concern  to  the 
experimenter,  and  since  the  concentration  varies  but  little  over  the 
cross-section  of  the  column. 

18*  l,ot  us  therefore  integrate  eq.  (15)  over  T  frfcm  T]_  to  Tg. 
Making  use  of  the  boundary  conditions  (17)  we  obtain  a  result  which 
can  be  written 


where  v  is  the  upward  transport  of  species  1,  given  by 


r  t,  i  M>c 

t  J  *  1 

o 


_  b  { 'A  N  n 

~Q  J  f  A.  I  VC,  -•  Jm  \  )  J  I  . 


Eq** (18)  states  the  obvious  fact,  that  r  is  constant  along  the  tube 
for  a  stationary  state.  *“  ~~— 

20*  Let  us  now  define  a  function 


A  M  -  Q  J  Aj-va/T 


We  see  that 


M? ;)=o, 


Hr)  =  , 


where  0/  is  the  total  transport,  given  by 


0  v:  k  I  pv  <rly 


h  ( ,J-  __ 

-:r  A  p  v  ci  T. 
Q  ~L  J 


Eq.  (19)  can  be  written 


■6  r* 


~  -q  J  AfJ)  gT 
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Integration  by  par.tfc,  with  tho  aid  of  (21),  gives 


~x  Tx 

__  <rc  IT  -<)  -  — .  I  F  •—  jt  -k.  f  Api>  </t 

7  -  °  d/  9/.  J  '  ■  (24) 

21*  If  one  now  goes  back  to  oq.  (15)  and  integrates  it  with 
respect  to  T,  one  obtains  an  expression  for  3c ,/ST  ,  which 

can  then  be  substituted  into  (24)*  Taking  into  account  the  boundary 
conditions  (17),  one  finds 

7 T  “  ~  *  f D  <J  ( V  5J  ~  ^  jp)  V  d  T  (25) 

and,  substituting  this  into  (24),  one  obtains. 

7  ~  ‘J{fT 

,  b  f7-  <rT'  J  T‘  (26) 

1  -v  r-  i  *  _  \  f 


M  dTJ  x^%dr  -  *  J  r- 


22,  If  one  noglects  the  small  variation  in  concentration  along 
a  section  zm  const.,  one  can  write  this  equation  as 

V  =  O-C,  +  P'  -,V  jp  ,  .27) 

where  H  ,  K,  and  N  are  positive  quantities  defined  by 


=  -  —  ( 
n  J 


nF*'r- 


K  *  K,, 

K<*  Hff'r, 

K  -  A.  (  >2- 
4  ~  QJ_'\yJ)4T, 
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and  the  concentrations  are  now  regarded  as  depending  only  on  z< 


23*  One  can  show  that  in  practical  cases  the  term  in  (27) 
involving  rT^c./ <y  will  generally  be  small  compared  to  the 

other  quantities  present.  Hence  we  shall  discard  this  term*  For 
convenience  we  shall  write  c  in  place  of  c  ,  1-c  in  place  of  eg* 
and  t  in  place  of  ~ct  ,  The  equation  tnen  becomes. 


C  r=.  a  e 


r  He  (»-<■)  ~  K  ~  . 


(a)  Equilibrium 


In  the  case  of  equilibrium  both  d  and  X  are  zero*  Eq*  (29)  then 
becomes 


H  c((-z)  —  o. 


If  we  lot 


This  can  bo  into grated  to  give 


=  e<  '>  u.  Const 


24*  Let  the  value  of  c  at  the  top  of  the  column  be  denoted  by 
,  that  at  the  bottom  by  c..  •  It  is  customary  to  define  the 

separation  factor  by 


C  <-%.  i./-  C.'l. 

J=  few'  <m> 

If  we  let  Se  denote  the  value  of  S  at  equilibrium,  then  it  follows 
from  (32)  that 

Ik  5e  =  *  L  }  Se  =  x°'L.  (34) 


Where  L  is  the  length  of  the  column*  An  important  feature  of  this 
result  is  that  the  value  of  the  equilibrium  separation  factor  is 
independent  of  the  concentration  at  the  bottom  of  the  column* 

25*  It  is  desirable  to  investigate  more  closely  the  dependence 
of  S0  on  the  variables  of  the  system*  In  terms  of  tho  integrals 
listed  in.  (28),  it  can  be  written 


lYtS  = 
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or,  by  uso  of  (2), 

ja ZkALkAAL  -  m 

«  'j'^T^-irDoT 

26#  If  ono  looks  for  the  value  of  a,  denoted  by  for  which 
lnSe  has  a  maximum  value,  one  finds  thr.t~it  satisfies  the  relation 


-L  ¥X>  Jr 

*.  / 

The  maximum  value  of  lnSe  is  then  given  by 


2h  = 


-/j  r  * 

-l  FetJr 


A*r 


This  depends  only  on  tho  wall  temperatures  T  and  T  #. 
is  rather  complicated,  however#  ^  ^ 


The  dependence 


27#  It  might  be  pointed  out  that,  on  the  basis  of  rough 
considerations  of  tho  behavior  of  the  various  quantities  in  (38), 
raising  T2  should  goncrally  increase  Sem#  However,  this  cannot  be 
done  indefinitely,  since  a  high  value  of  &T  will  lead  to  parasitic, 
or  local  convection,  resulting  in  mixing  of  the  liquid  and  lowering 
of  the  separation  factor.  It  appears  that  the  choice  of  optimum  wall 
temperatures  can  best  be  made  on  the  basis  of  experiment,  rather  than 
theory# 

28#  Returning  to  (36),  we  see,,  that,  for  fixed  wall  temperatures 
Tj  and  Tg,  the  dependence  of  SQ  on  a  oan  be  expressed  in  the  form 

jj=  m 

c  /  a 

where  kliand  ^  are  eonstants  depending  on  the  temperature  •  A  similar 
relation  was  given  by  YJaldmonn*-  for  gases# 

^  29#In  plate  l,the  ^crosses  indicate  experimental  values  of  InS 
for  a  number  of  values  of  a,  as  determined  by  Dr.  Philip  H.  Abelson 
for  a  particular  set  of  conditions#  The  curve  represents  In  S0  as 
givon  by  (39),  k^  and  kg  having  been  chosen  to  fit  the  experimental 
data#  It  will  be  seen  that  the  agreement  is  good.  However,  the 
agreement  may  be  fortuitous  to  some  extant  for,  as  Dr#  Abelson 
pointed  cut,  in  the  wxperiments  what  was  kept  constant  was  not  the 
temperature  of  tho  wall  in  contact  with  tho  diffusing  liquid,  but 
rather  the  temperature  of  the  heating  or  cooling  substance  in  contact 
with  the  other  sido  of 
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the  wall.  Vjhen  the  distance  a  is  varied,  the  temperature  of  the 
diffusing  liquid  at  the  wall  will  vary  somewhat  in  spite  of  the  fact 
that  the  heating  and  cooling  agents  are  kept  at  fixed  temperatures. 

30.  If  we  wish  to  have  a  solution  for  the  concentration  c  at 
equilibrium,  we  can  obtain  it  in  a  convenient  form  by  expanding  it  in 
a  power  series  in  z.  Substituting  into  eq.  (30)  and  equating  to 
zero  the  coefficients  of  each  power,  we  get 

=  e0  -f-  r  (i-  -'o)c^2  ^2^  •••)  >  (40) 

provided  we  take  z  =  0  at  the  bottom  of  the  column. 

31.  If  the  solution  is  desired  in  a  closed  form  rather  than  as  a 
power  series,  it  can  be  obtained  from  (32).  The  result  is 


C  ( z )  r  C  £  u 

"  l^Ce 

(b)  Stationary  flow 


(40a) 


32.  Vie  next  consider  the  case  in  which  there  is  a  steady  upward 
flow  of  liquid  through  the  column  at  such  a  rate  that  the  concentration 
at  every  point  remains  constant  with  time,  a  condition  which  might  be 
referred  to  as  "stationary".  In  this  case  eq.  (29)  applies,  withO 
and  r  both  constant, being  the  mass- of  liquid  flovdng  upward  per 
unit  time  through  any  cross-section  of  the  column,  while  f  is  the  mass 
of  the  desired  constituent  crossing  any  section  per  unit  time. 

These  two  quantities  are  related  to  each  other  by  the  equation 

7U,  (41) 

where  %  is  the  value  of  c  at  the  point  where  the  enriched  material  is 
withdrawn.  If  one  is  operating  a  single  column  this  \vill  be  the  same 
as  C^.  i  but  if  the  column  is  part  o-f  a  pyramid,  will  be  the 
concentration  at  the  top  of  the  pyramid f  Let  us  for  the  present  impose 
no  restriction  on  c  • 

Eq.  (29)  can  now  be  written 

, 'v' cm  2C~°  •+■  Hc(l-c)  -  K  do,  (42) 


If  we  define  as  in  (31)  and  let 

v  :  t 

H  , 

the  equation  becomes 
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This  can  be  solved  most  conveniently  by  means  of  a  power  series  in  z . 
In  this  way  one  obtains  as  the  solution 


<5  =  «n  + 


j^Co(l-0o)  -  T  (®m“C0)||jy^.  4"^  (^”2co  "f*  'v)^-2 


putting  z  s  L  gives  us  C  ,  and  f  rom  this  #e  get  the  separation  factor 


s  s  l^if,  '  -AfW-Oe  f\  ,  [*,.  j /in,- 

r  by^—rz — J  ; (46) 

33.  If  we  are  dealing  with  the  case  in  which  cm  -  c  ,  r  we  can 
put  the  result  into  a  more  convenient  form.  At ‘the  top  pf  the  column 
(45)  then  becomes 

c^_=  Co-f-f  c0(l-co)  -v(C-f.”  ojfa  L  ^§(l-2c0  • 


y  (h  ir!^  c<-~ 


/-  e, 


we  solve  this  equation  for  c,  we  obtain 

+ 


Vj 

J.,  -  • 

(46) 

,  we 

can 

the 

column 

,  2  l 

%'  -f.  • 

j  ■ 

(47) 

(48) 

(49) 

c_|_-  c0+co(l-co)  CX  L  -f  (l-2c0-y)<^2  L2]  -f-»  *  •  (48) 

and  from  this  vre  find 

S  r  1  +  t\L  f  |  (1-V1)^  21_^V  4‘,*‘  (49) 

34.  It  will  be  seen  by  a  comparison  of  (49)  with  (34)  that  when 
liquid  is  xui thdrawn  continuously  from  the  top  of  the  column  the 
separation  factor  is  less  than  that  at  equilibrium;  the  difference 
will  be  very  small  however  if  Se  is  nearly  equal  to  unity  QjL  small). 

APPROACH  TO  EQUILIBRIUM 

35.  If  one  is  interested  in  following  the  process  by  which  the 
column  goes  from  its  initial  state  to  the  condition  of  equilibrium, 
one  must  replace  eq.  (29)  by  the  more  general  equation  involving 

the  time.  This  equation  can  be  obtained  from  eq.  (12)  in  much  the  same 
way  as  (29)  was  obtained  from  (13).  Such  an  equation  was  derived  by 
Bardeen.'  If  we  again  omit  terms  which  are  small  in  practical  cases, 
the  equation  can  be  written  in  the  .form 


y  ” k  As.  -  H  3l-I 

-3  t  5  3  ' 

If  we  define  c^as  in  (31)  and  let 


_.[  c  (■/-<;)]  -tK  cij 


this  becomes 


j  -  Ai  j  -  yo.  b 

b\  y  H 


0  -11-  -* 
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36.  To  solve  this  equation,  we  look  for  a  solution  which  is  a 
power  series  in  z,  with  coefficients  which  are  functions  of  t. 

Let  us  set 

c  z  30  (t)  +  ifi  ^  -f  •••  (53) 

Substituting  this  into  (52)  and  equating  coefficients  of  corresponding 
powers  of  z  on  both  sides  of  the  equation,  we  obtain  a  series  of  relat- 


(j  *  h)  +  y  \ 

r -  - -*• /.*>. (/-<* ^ - yf] + 6 y3  > 

7  =  ~  Y?  /^Yv. 

rime  denotes  differentiation  with  respect  to  time. 


>  (54) 


where  a  prime ' denotes  differentiation  with  respect  to  time, 
be  seen  that  two  of  the  functions,  X  0  and  ^  can  be  taken 
■arbitrarily.  The  others  will  then  be  determined  by  the  equa 
terms  of  these  two.  One  finds,  for  example, 

■The  upward  transport  X,  given  by 


It  will 


quations  in 


E  -  He  (/-«)-  K  '|-X,  ? 


is  found  to  be 

r  -  M  \  Y  (/-  y  Vy  .  J/ h  Yt  x  y  v,.\,-±s3  Y 
n  jjb  {/  V0j  cr/  >  of 4  /  4  df  Y~/3  j 4  # -  •  j* 

37.  The  functions  X  and  C)  q  will  be  chosen  in  any  perticular 
problem  so  as  to  satisfy  the  initial  conditions  and  the  boundary 
conditions  at  the  top  and  bottom  of  the  column. 

38.  This  method  of  solution  is  much  simpler  than  the  usual  method 
involving  a  Fourier  series,  and  it  should  work  satisfactorily  for  the 
cases  encountered  in  practice.  An  example  of  its  application  will 
next  be  given. 

(a)  Fixed  concentration  at  bottom  of  column. 

39.  he  consider  the  case  of  a  column  of  the  type  previously 
discussed  which  is  operated  by  having  its  lower  end  connected  to  a 
large  reservoir  or  to  a  stripper  column  so  that  the  initial  concen¬ 
tration  is  maintained  there  at  all  times.  i!e  have  then,  by  taking 

z  a  0, 


5  c  ( t)  — ; 


t  •-  r  -  >  \  t 


-  C, 
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whence  it  follows  from  (54)  that 

K?.  '  ^  l  f  ^  ”  0  c  ) 

Wd'-^coytfl+KlK’ 


$¥=  yrlWt-  .  J 

)•  In  this  case,  eq.  (57)  gives  for  r  ' 

r  =  H  [co  u-<c)-yr*  faty-xii 


At  the  top  of  the  column  (z  ±  L)  we  have 

+ b oO c-)~ %] '/>i  fc  ItU,  LX jf /  bf y/ (%r ‘ 


Hi 


or,  s ince 


M-.:  t>  L 


we  can  write 

/* f^r/s.-^>vLJ+"  <6! 

41.  On  the  other  hand,  at  the  top  of  the  column,  we  have  the 
rate  of  change  of  concentration  related  to  the  transport, 

r+  -  Mn  4±+  (6- 

+  n  </t 

where  JL  is  the  mass  of  the  liquid  in  the  reservoir,  while  c+is 
given  by 


v-i.'  * 


Substituting  (65)  and  (63)  into  (64)  we  get  an  equation  for  )()  ; 

yj-  4  L  S/  f \!j 

Neglecting  higher  order  terms,  one  integrates  this  to  get 

V-c.0-c„)[uc-^]  (67) 

where  the  relaxation  time  ©  is  given  by 

6  4iH,  [>  %  -  -  ok l]h%  fnm -<•»)  ulV  (68) 
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to  the  approximation  considered 


42.  The  solution  for  the  concentration  can  be  written 


C^o  + 


-  u  0 -  C o  )  (l  -  "f  +  ( hr  t  hr  J  "  j  ~h.  .  . 


From  this  one  finds  for  the  separation  factor 


5  -  /  r(  /  -  e  -Vs] )[<A  L-t  <>/x- cv  e  'V*jrL-sj+. . . 
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